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A  Simple  Model  for  the  Initial  Phase  of  a  Water  Plasma  Cloud  about  a  Large 

Structure  in  Space. 

D.E.  Hastings,  N.  A.  Gatsonis  and  T.  Mogstad 
Dept  of  Aeronautics  and  Astronautics, MIT,Cambridge,  MA  02139 

Abstract 

Large  structures  in  the  ionosphere  will  outgas  or  eject  neutral  water  and  perturb  the  ambient 
neutral  environment.  This  water  can  undergo  charge  exchange  with  the  ambient  oxygen  ions  and 
form  a  water  plasma  cloud.  Additionally,  water  dumps  or  thruster  firings  can  create  a  water  plasma 
cloud.  A  simple  model  for  the  evolution  of  a  water  plasma  cloud  about  a  large  space  structure  is 
obtained.  It  is  shown  that  if  the  electron  density  around  a  large  space  structure  is  substantially 
enhanced  above  the  ambient  density  then  the  plasma  cloud  will  move  away  from  the  structure.  As 
the  cloud  moves  away  it  will  become  unstable  and  will  eventually  break  up  into  filaments.  A  true 
steady  state  will  exist  only  if  the  total  electron  density  is  unperturbed  from  the  ambient  density 
When  the  water  density  is  taken  to  be  consistent  with  Shuttle  based  observations  the  cloud  is  found 
to  slowly  drift  away  on  a  timescale  of  many  tens  of  milliseconds.  This  time  is  consistent  with  the 
Shuttle  observations. 

1  Introduction 

Since  the  dawn  of  Shuttle  flights  to  the  ionosphere  it  has  been  realised  that  the  ambient  iono¬ 
sphere  is  strongly  perturbed  by  the  presence  of  a  large  body  in  spacejGreen  et  al . ,  1985] .  Several 
experiments  on  board  the  Space  Shuttle  have  been  devoted  to  the  measurement  of  the  plasma 
parameters  in  the  vicinity  of  the  spacecraft[Murphy  et  al. ,  1986;  Reasoner  et  al.,  1986).  It  has  been 
found  that  the  Shuttle  or  any  body  in  space  affects  the  ambient  environment  through  the  inter¬ 
action  of  its  associated  contaminant  cloud  with  the  environment.  The  contaminant  cloud  arises 
from  surface  outgassing,  leaks  from  life  support  systems,  thruster  firings  and  water  dumps.  The 
contaminant  cloud  can  undergo  several  different  types  of  interactions  with  the  ambient  environ¬ 
ment  ranging  from  physical  interactions  such  as  momentum  transfer  to  chemical  interactions  such 
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as  charge  exchange.  The  result  is  that  the  self-consistent  environment  around  a  large  body  in 
space  may  be  significantly  different  from  the  ambient  environment.  We  choose  to  put  the  empha¬ 
sis  on  large  bodies  since  large  bodies  can  have  an  associated  contaminant  cloud  which  makes  an 
0(1)  change  to  the  ambient  environment.  By  an  0(1)  change  we  mean  that  the  composition  of 
the  ion  density  is  substantially  changed.  For  example,  observations  around  the  ShuttlejPickett  et 
al.,1985]  indicate  that  the  observed  environment  differs  considerably  from  the  expected  ambient 
environment. 

It  is  important  to  study  the  self-consistent  environment  around  large  space  vehicles  such  as  the 
Shuttle  for  several  reasons.  The  major  reason  is  to  understand  the  noise  that  will  be  observed  on 
any  sensors  carried  on  the  vehicle[Pickett  et  al.,1985].  This  will  directly  affect  the  utility  of  the 
sensors  for  purposes  of  observation.  Another  important  reason  is  to  build  understanding  of  the 
space  plasma  environment  and  some  of  the  processes  that  can  occur  in  it.  Finally  study  of  the 
self-consistent  environment  will  help  in  understanding  such  basic  design  issues  as  how  to  mitigate 
spacecraft  charging  in  low  earth  orbit. 

In  this  paper  we  develop  a  simple  model  for  the  motion  of  a  water  plasma  cloud  around  a 
large  structure  in  the  ionosphere.  In  section  2  we  derive  the  equations  that  describe  the  large  scale 
plasma  flow  around  the  body.  In  section  3  we  obtain  a  analytic  solution  for  the  case  where  there 
is  no  electrostatic  modification  of  the  motional  potential.  This  case  corresponds  to  a  plasma  cloud 
which  does  not  enhance  the  ambient  density.  In  section  4  we  perform  a  multiple  scale  analysis  and 
obtain  the  equations  for  the  field  line  averaged  densities  and  potential.  For  the  case  where  the  cloud 
density  exceeds  the  ambient  density  we  show  that  the  equations  of  motion  do  not  allow  steady  state 
solutions  and  suggest  that  the  cloud  may  be  unstable  to  the  growth  of  E  x  B  gradient  instabilities 
We  then  solve  the  model  equations  numerically  with  an  Flux  Conserving  Transport(FCT)  numerical 
method  and  show  the  long  time  behaviour  of  the  plasma  cloud.  We  obtain  the  ion  residence  time 
for  a  range  of  cloud  densities.  For  densities  typical  of  what  is  observed  around  the  Shuttle,  the 
residence  times  are  consistent  with  the  measurements.  Finally  in  section  5  we  conclude  with  a 
discussion  of  the  significance  of  these  results  for  measurements  from  large  space  vehicles  such  as 
the  Shuttle.  Although  the  equations  of  the  water  plasma  cloud  are  generally  applicable  we  choose 


to  concentrate  on  the  dynamic  behavior  of  an  initial  water  plasma  "puff”  which  could  arise  from 
an  MPD  thruster  firing  or  a  plasma  source. 

2  Derivation  of  the  equations  for  the  plasma  cloud 

Plasma  clouds  in  the  ionosphere  have  been  studied  both  theoretically  and  experimentally  for 
a  number  of  years(Perkins  et  al.,1973;  Zabusky  et  al.;  1973;  Mitchell  et  al.,  1985].  The  major 
motivation  has  been  to  understand  the  dynamics  of  such  clouds  since  they  have  been  used  for 
tracing  the  magnetic  field  lines  and  they  may  interfere  with  communications.  In  this  work  we 
shall  follow  a  similiar  path  to  obtain  the  equations  for  the  plasma  and  self-consistent  electrostatic 
potential  around  a  large  body. 

We  work  in  a  frame  attached  to  the  space  structure  and  start  by  assuming  that  the  space 
structure  is  emitting  neutral  water  into  an  ambient  background  of  oxygen  ions  and  neutrals.  The 
spacecraft  is  taken  to  be  moving  at  orbital  velocity  with  respect  to  the  ambient  environment^  8 
km/sec  for  low  earth  orbit).  For  an  ambient  environment  whose  neutral  density  is  of  the  order  of 
108  cm-3  the  mean  free  path  of  the  water  molecules  is  many  kilometers.  Hence  for  length  scales 
of  up  to  a  few  kilometers  we  can  take  the  water  molecules  as  freely  expanding  with  their  thermal 
velocity  and  write  the  steady  state  neutral  water  density  as 

=  n°HjO{r0/r)2  (1) 

This  assumption  about  the  neutral  density  is  reasonable  for  low  water  densities  (water  density  much 
less  than  the  ambient  neutral  density)  but  will  have  to  be  modified  for  higher  densities.  However 
since  we  wish  to  concentrate  on  the  plasma  dynamics  we  shall  leave  more  detailed  modelling  of 
the  neutrals  for  future  work.  In  Eq.  (1)  we  have  approximated  the  structure  as  a  sphere  of  radius 
r0  and  taken  r  as  the  radial  distance  from  the  center  of  the  sphere.  This  approximation  will  be 
reasonable  for  any  body  shape  as  long  as  we  are  several  body  dimensions  from  the  body.  So 
for  example,  for  the  Space  Shuttle  whose  longest  dimension  is  of  the  order  of  40  meters  can  be 
reasonably  approximated  in  this  manner  at  distances  of  the  order  of  100  meters  and  beyond.  The 
density  is  the  density  of  water  at  the  surface.  The  assumption  that  the  neutral  species  is 
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water  is  based  on  observations  from  the  Space  Shuttle  where  the  dominant  contaminant  species 
was  observed  to  be  HjO+ [Caledonia  et  al.,1987]. 

The  ionic  water  may  be  formed  by  the  charge  exchange  reaction 

0  +  +  HiO  O  +  H20+  (2) 

where  in  the  frame  moving  with  the  space  structure  the  oxygen  ions  sweep  through  with  an  energy 
of  up  to  5  eV.  At  this  energy  the  reaction  rate  kcz  has  been  measured  and  is  kcx  —  6  x  10-1C 
cm3/sec[Murad  and  Lai,  1986a]  and  for  0+  the  charge  exchange  collision  frequency  is  i>cz  — 

We  assume  for  simplicity  that  the  body  is  moving  perpendicular  to  the  magnetic  field  with 
velocity  in  a  fixed  frame  of  V .  We  define  a  wind  velocity  in  the  moving  frame  as  =  -V  and 
then  define  a  cartesian  coordinate  system  (z,y,  z)  in  the  moving  frame  by  taking  the  distance  x  to 
be  along  the  direction  of  -VWt  the  distance  y  to  be  along  the  direction  of  -V  x  B  and  the  distance 
z  to  be  along  the  direction  of  the  magnetic  field  B.  We  note  that  the  negative  y  direction  is  along 
the  direction  of  the  unshielded  motional  electric  field  that  will  be  seen  in  the  far  field  from  the 
moving  frame.  The  frame  is  illustrated  in  Figure  1 

We  take  the  oxygen  and  water  ions  to  be  mainly  moving  across  the  magnetic  field  due  to  their 
large  gyroradii  while  the  electrons  are  taken  to  dominantly  flow  along  the  magnetic  field  due  to 
their  small  gyroradii.  This  assumption  on  the  ion  motion  is  true  on  the  large  scale  where  the  ions 
have  a  chance  to  complete  their  gyroorbits.  In  Caledonia  et  ai  1987]  the  ions  were  taken  to  be 
almostly  perfectly  shielded  from  the  motional  electric  field  and  hence  in  the  moving  frame  the  only 
motion  available  to  them  was  parallel  motion  In  this  work  the  potential  is  being  calculated  in  a 
self-consistent  manner  so  that  the  ions  can  move  accross  the  magnetic  field.  The  perpendicular 
velocity  of  the  oxygen  ions  can  be  obtained  from  the  steady  state  momentum  balance  equation  Krall 
and  Trivelpiece,1973j  and  is 

Cot  -  (kQ4 -Fl  -  * X  6)  (1  -  K2,.  )  ; 

where  6  is  the  unit  vector  in  the  magnetic  field  direction  with  k,,.  being  the  Hall  parameter  given 


*0+  —  ^O^/i^O+.HiO  +  Vex).  (-1) 

and  Q0+  being  the  oxygen  cyclotron  frequency  while  t'o+,//Jo  ‘s  the  oxygen  ion-neutral  water 
momentum  exchange  collision  frequency.  The  vector  F±  is  given  by 

F±  =  ~{cTi/eB)V x  lnn0+  4-  ( ce/B)E  +  i/0+ x/Ho+  ■  (5) 

In  Eq.  (5)  T,  is  the  ion  temperature  which  we  take  to  be  constant  and  the  first  term  is  the  force 
due  to  the  ion  pressure,  E  is  the  electric  field  seen  from  the  moving  frame  and  gives  the  electric 
force  on  the  ions  in  the  moving  frame,  while  vh-iO  in  the  third  term  is  the  velocity  of  the  water 
and  this  term  expresses  the  force  on  the  oxygen  ions  due  to  collisions  with  the  water  neutrals.  In 
the  derivation  of  the  force  on  the  oxygen  ions  the  dominant  collisions  for  these  ions  has  been  taken 
to  be  collisions  with  the  neutral  water  molecules(rate  constants  1.2  x  10'9  cm3/sec).  Collisions 
with  oxygen  neutrals(rate  constants  3  x  10-9  cm9/sec)  have  been  neglected  for  simplicity  since 
we  anticipate  little  differential  velocity  between  oxygen  neutrals  and  ions.  In  Eq.  (3)  the  first  term 
gives  rise  to  the  well  known  Pedersen  conductivity  while  the  second  contains  the  E  x  B  drift  and 
the  Hall  current.  We  note  that  the  Hall  parameter  takes  very  large  values  for  parameters  typical 
of  the  ionosphere(i/0+  H^0  ~  2  x  10~3  sec-1  and  Q0+  —  1-8  x  I02  sec-1  so  that  /c0+  ~  9  x  10*). 
This  suggests  that  the  dominant  motion  for  the  oxygen  ions  will  be  the  E  x  B  drift.  Hence  we  can 
write  the  continuity  equation  for  the  oxygen  ions  as 

^  x  6  - Vn0+  =  -n0+i/CI  (6) 

where  the  electric  field  in  the  moving  frame  is  E  —  -V<£  and  the  boundary  condition  on  Eq.  (6)  is 
that  for  jr!  — *  oo  the  oxygen  density  approach  the  ambient  density. 

The  electrons  will  mainly  flow  along  the  magnetic  field  but  will  also  have  an  E  x  B  drift  across 
the  field.  In  addition  electrons  will  be  lost  as  a  result  of  dissociative  recombination  of  the  water 
ion.  This  occurs  with  a  rate  krecom  which  we  shall  specify  later.  Hence  we  can  write  the  electron 


continuity  equation  as 
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where  IY.  is  the  parallel  electron  flux  given  by 


and  the  potential  tpe  is  defined  as 
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In  Eq.  (9),  <f>  is  the  self-consistent  potential  which  is  modified  by  the  electron  pressure.  The  parallel 
current  is  related  to  the  parallel  electric  field  through  the  resistivity (r?«.) 


~eV‘«  =  rT  (_VHV'e) 


where  the  resistivity  is  given  by 
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In  Eq.  (8)  the  collision  frequency  ut  is  the  electron  momentum  exchange  frequency  due  to  collisions 
with  the  neutrals  or  the  ions.  The  physical  content  of  this  equation  is  that  the  electrons  flow  along 
the  field  lines  and  are  swept  across  it  by  the  self  consistent  potential  in  the  plasma  cloud  We 
note  that  for  any  steady  state  to  exist  to  this  equation  it  is  necessary  that  there  be  a  convection 
of  electrons  across  the  field  lines.  This  implies  that  the  potential  contours  must  be  open.  This  is 
because  closed  potential  contours  will  cause  the  plasma  to  rotateiDaly  and  Whalen,  19791  but  will 
not  lead  to  any  net  loss.  This  implies  that  the  symmetry  breaking  part  of  the  electrostatic  potential 
will  be  crucial  in  determining  the  motion  of  electrons.  This  part  of  the  potential  will  arise  from 
the  fact  that  the  system  is  moving  and  hence  sees  a  motional  potential 

Once  we  have  the  equations  for  the  oxygen  ions  and  the  electrons  the  equation  for  the  water 
ions  follows  from  the  requirement  that  the  plasma  remain  quasineutral, 


nH-)0+  -  ne  ~  nO»  (  17 

As  an  alternative  to  using  Eqs.  (6),  (7),  and  (12)  we  could  use  Eqs.  (6)  and  (12)  coupled  with  th 
equation  for  the  water  ions 


d"H2o+  _  c_ 
dt  B 


V ^<p  x  b  -  Vrij/j  —  n0»i/cx  nenHi0 
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Finally  we  need  an  equation  to  determine  the  self-consistent  potential  .  This  equation  comes 
from  the  requirement  of  quasineutrality  for  all  times.  This  means  that  in  addition  to  Eq.  (12)  we 
require  that  we  have  charge  conservation (V  •  j  =  0  where  j  is  the  current  flowing  in  the  plasma) 
We  take  the  perpendicular  current  to  be  given  by  the  perpendicular  ion  flows(  Eq.  (3)  and  the 
analogous  equation  for  the  perpendicular  water  velocity)  and  the  parallel  current  to  be  given  Dy 
the  electron  flows(  Eq.  (10))  and  obtain  the  equation[Drake  and  Huba,  1986] 
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In  this  equation  for  charge  conservation  the  perpendicular  ion  diffusion  coefficient  is 
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where  the  Hall  parameter  for  the  water  ions  is 


kH20+  -  ^HjO  +  Kv/ijO  +  .O  ~  vcj)  (lb) 

We  note  that  in  this  expression  all  the  E  x  B  drifts  have  cancelled  out  since  they  do  not  give  rise  to 
any  net  current.  The  first  term  on  the  RHS  of  Eq.  (14)  arises  from  the  Pedersen  current  due  to  the 
electric  field,  the  second  term  due  to  the  diffusion  of  ions  down  the  density  gradients,  the  third  due 
to  the  Hall  current  while  the  fourth  is  the  electron  flow  to  balance  the  ions.  This  equation  must  be 
solved  subject  to  the  boundary  condition  that  in  the  far  field  the  electric  field  in  the  moving  frame 
is  the  motional  electric  field.  This  requires  that  for  |f]  — *  oo  we  have  =  -V  x  B/c.  We  can 
interpret  the  first  two  terms  in  the  potential  equation(  Eq.  (14)  as  giving  the  shielding  from  the 
plasma  cloud  while  the  last  term  is  the  charge  neutralization  by  the  parallel  flow  of  electrons  from 
far  away  in  the  ionosphere. 


3  Analytic  Solutions  of  the  Cloud  Equations 


The  set  of  equations  for  the  potential  and  densities  cannot  in  general  be  solved  analytically 
However  it  is  possible  to  obtain  a  special  solution  and  draw  some  general  conclusions  from  the 
structure  of  the  equations. 

If  the  density  of  the  cloud  is  very  low  then  we  expect  the.t  the  motional  potential  will  be  oui\ 
weakly  shielded  by  the  cloud.  In  this  case  we  approximate  the  electric  field  everywhere  L\  •:.< 
boundary  condition  on  Eq.  (14)  and  write  -  -V  x  B/r  for  ail  space  if  we  require  -py 

4>  =  0(rl)  (i.  e.  we  do  not  allow  any  parallel  electric  fields)  then  iri  the  absence  of  rer  rno;  n  v. . 

we  have  from  Eq.  (7)  in  steady  state  that 

-  V'  •  V  =  0.  17 


which  has  only  two  possible  solutions,  either  that  nt  -  njmiltru  everywhere  or  that  n,  varies  r 
only  one  direction  in  space  which  is  physically  unreasonable  The  oxygen  equation  is 
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which  has  the  solution 
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This  shows  the  product  of  the  water  column  density  traversed  by  a  O  '  i..n  and  the  ,  h  tri-  -  \  ...v.i 
cross  section.  Hence  the  electron  density  is  constant  everwher-  while  •.  yg>  :  \  v;  .■  • . 

and  partially  replaced  by  the  water  density.  The  contours  of  .r.t  a  c  i  . 

to  the  ambient  electron  density  are  shown  in  Figure  J  f. < r  a  •  .  ,  ..  ii  ,  .v-e  ;  a:,  •  .utga.ss.ng  • 
which  is  emitting  water  from  a  central  region  of  ICO  meter ••  r  i.ia-  U.-  .  b  ,<•.>«»  •  ,  K  • ; . .  - . 
large  length  scales  since  there  is  evidence  from  PDF  m<  .isar*  i:  i  •  •  SI,  nr  >  .  f  w.y  • 

ions  out  to  a  kilometer  from  the  Shuttle  G  Murphy,  private  .  .y.,.n,  i  We  see  ••  a 

most  of  the  depletion  for  the  oxygen  density  occurs  on  a  lr.<  ,;,1»  ;  • y  .re  a;.d  f,  r  * 

case  somewhat  less  than  half  of  the  oxygen  ions  are  convert e«i  t..  water  •.■ns  H'-nre  Mr  this  y-. 
the  water  plasma  cloud  is  actually  a  depletion  relative  to  the  ammei.t  . > i t ■  <  tron  density 


4  Derivation  and  Numerical  Solution  of  the  Averaged  Cloud 


Equations 

Since  the  plasma  flow  along  the  magnetic  field  is  extremely  rapid  compared  to  the  perpendicular 
flow  we  can  simplify  the  equations  describing  the  plasma  by  doing  a  formal  multiple  scale  analysis. 
Physically,  we  expect  the  cloud  to  take  a  cigar  shape  along  the  magnetic  field.  The  equations  for 
the  field  line  averaged  density  and  field  line  averaged  potential  can  be  obtained  by  averaging  the 
density  and  potential  equations  where  the  averaging  operator  is  given  by 


F  =  — ^ —  f  f  — 

f  dl/B  J  B 


/  dlj  B  J  J  B 

Hence  from  this  point,  Nt,  Nh2o+  .  No  +  >  $  refer  to  the  field  line  averaged  electron  density, water  ion 
density  etc.  The  field  line  average  operator  annihilates  the  parallel  terms  in  the  electron  density 
and  potential  equations  ard  gives  for  the  density  equations  (from  Eqs.  (6)  and  (7)) 

^-§  V±*xi.VWo,=  -AT0^,  (20) 
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while  the  potential  equation  (  Eq.  ( 14) )is 
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In  obtaining  these  equations  we  have  used  the  multiple  scale  assumption  to  say  that  the  densiti?s 
and  potential  have  a  strong  dependence  on  the  perpendicular  direction  and  a  weak  dependence  on 
the  parallel  direction.  This  is  a  reasonable  assumption  for  the  plasma  quantities  but  cannot  be 
applied  to  the  neutrals.  Therefore  in  Equations  (20),  (21),  and  (22)  we  have  evaluated  the  neutral 
water  density  in  the  plane  2  =  0.  This  assumption  enables  us  to  obtain  a  closed  set  of  equations 
for  the  field  line  averaged  quantities  and  still  contains  the  appropriate  physics. 


This  set  of  equations  describes  the  dynamic  evolution  of  a  water  plasma  cloud  as  seen  from  a 
moving  frame  and  with  the  potential  calculated  self-consistently.  We  shall  use  these  equations  to 
study  the  evolution  a  water  plasma  cloud  where  the  initial  density  may  be  larger  than  the  ambient 
electron  density.  As  pointed  out  in  Caledonia  et  al.[1987]  such  an  enhanced  density  cannot  arise 
self-consistently  solely  by  chemical  processes  such  as  charge  exchange.  Therefore  in  those  cases 
where  the  water  plasma  density  exceeds  the  ambient  density  we  assume  that  the  water  plasma  has 
been  injected  at  the  initial  time  due  to  some  mechanical  mechanism  such  as  the  firing  of  a  plasma 
source. 

This  set  of  field  line  averaged  equations  for  a  plasma  cloud  around  an  space  structure  still  cannot 
be  solved  analytically  but  it  can  be  shown  that  the  only  steady  state  solution  allowed  is  the  one 
previously  discussed  where  the  electron  density  is  uniform  everywhere.  The  proof  of  this  statement 
is  given  in  the  appendix  and  follows  a  similiar  proof  in  work  by  Dungeyj  19581.  Since  in  the  far 
fhld  the  electron  density  must  match  the  ambient  density  we  conclude  that  plasma  clouds  whose 
density  exceeds  the  ambient  density  cannot  be  steady  state  structures.  Therefore  the  enhanced 
densities  seen  around  the  Shuttle[Pickett  et  al. ,19851  must  be  transient  phenomena.  The  fact  '.hat 
a  plasma  cloud  whose  density  exceeds  the  ambient  will  continously  evolve  is  also  suggested  by 
the  following  physical  arguments:  such  a  plasma  cloud  will  have  a  density  distribution  which  will 
probably  peak  towards  the  center  of  the  cloud.  The  moving  cloud  will  have  Pedersen  currents  in 
it  which  act  to  modify  the  imposed  motional  field.  These  Pedersen  currents  are  density  dependent 
since  they  arise  from  the  bulk  motion  of  ions  due  to  the  potential  gradients  Therefore  we  expect 
and  Eq.  (22)confirms  this,  that  the  highest  density  regions  cf  the  cloud  will  modify  the  motional 
potential  most  effectively.  Hence  the  low  density  cloud  edges  will  see  a  self  consistent  E  *  B  drift 
different  than  the  higher  density  regions.  This  leads  to  a  velocity  shear  through  the  cloud  which 
causes  it  to  continously  distort  Therefore  no  steady  state  will  exist.  The  timescale  on  which  this 
distortion  will  take  place  will  be  a  function  of  density  It  ran  be  estimated  on  the  basis  of  the 
arguments  given  by  Perkins  et  al.  1973;  if  we  assume  that  the  dominant  effect  of  the  cloud  is  to 
shield  the  potential.  If  the  cloud  has  density  nr  »  n.imb,ent  then  the  electric  field  in  the  cloud  is 
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Ev  =:  namW<nl/ne£b  where  Eq  =  V  x  B/c  =  £bv  From  V  x  £  =  0  we  obtain 

Ex  —  —  y(t*amA««nl /tlc  )  (dnc/dx)  Bo- 
The  cloud  velocity  is  =  c£  x  B/ B2  which  gives 

V*  —  —  V  (flambient  /  He) 


and 

vv  ~  -V(ynamiUrxt/n2e){dne/dx). 

We  see  that  in  the  moving  frame  the  cloud  will  drift  backwards  with  a  velocity  which  depends 
on  the  density  and  will  drift  sideways  with  a  velocity  which  depends  both  on  the  density  and  on 
how  sharp  the  density  gradient  is  in  the  direction  of  the  neutral  wind.  From  these  estimates  we 
can  see  that  distortion  will  occur  on  a  timescale  of  Ln/(^(nami«nt/”c))  where  Ln  is  the  density 
length  scale.  For  the  case  of  the  Space  Shuttle  if  we  take  the  density  length  scale  to  be  the  Shuttle 
length(~  50  meters)  and  take  a  water  ion  cloud  of  107  cm-3  in  an  ambient  density  of  105  cm  3 
then  for  a  velocity  of  8  km  per  second  the  distortion  time  is  0.6  seconds. 

We  have  solved  the  averaged  cloud  equations(Eqs.  (20),  (21),  and  (22))  for  a  5km  by  5km 
square  region  using  a  uniform  101  by  101  cartesian  mesh.  The  center  of  the  square  was  fixed  on 
the  structure. With  «//3o  +  (z >!/)  and  no+(z>y)  g>ven  at  some  time  we  calculate  the  potential  from 
Eq.  (22)  and  use  that  result  to  obtain  the  flow  field  velocities.  Then  that  information  is  used  in 
Eqs.  (20)  and  (21)  to  advance  the  densities  by  one  timstep. 

Successive  Point  Over  Relaxation(SPOR)  was  used  to  solve  the  elliptic  current  balance  equation 
This  scheme  introduces  a  relaxation  parameter  into  the  equation  and  which  then  allows  an  iterative 
solution  of  the  problem[Roach,1976] .  The  derivatives  were  approximated  by  second  order  finite 
differences  and  iterations  were  carried  out  until  the  maximum  error  was  less  than  10  ~5  between 
twosucessive  iterations.  Neumann  boundary  conditions  were  applied  in  both  the  x  and  y  directions 
We  required  the  potential  to  match  the  following  boundary  conditions  in  the  far  field 
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in  the  V  x  B  direction  and 

Vj>  =  0 

in  the  direction  of  motion.  The  potential  of  the  central  structure  waa  taken  to  be  4  Volts  which  is 
consistent  with  measurements  taken  from  the  Shuttle. 

For  the  convective  equations  we  used  the  two-dimensional  flux  correction  method(FCT)  of 
Zalesak[l979],  The  high  order  scheme  was  a  leapfrog-trapezoidal[Roach,1976]  with  fluxes  calculated 
with  the  flux  formulae  developed  by  Zalesak[1984j  while  a  donor-cell  scheme  [Roach,  1976]  was  used 
for  the  lower  scheme  to  complete  the  FCT  algorithm.  Symmetric  boundary  conditions  were  applied 
on  the  density  distributions  and  the  flux  limiter  was  applied  on  every  iteration.  No  assumptions 
were  made  about  the  symmetry  of  the  solutions  so  that  the  computations  were  carried  out  through 
the  whole  mesh.  The  oxygen  density  approached  the  ambient  ion  density(  taken  to  be  2  x  10s 
cm'5)  at  the  edges  of  the  square.  The  water  ion  density  was  assumed  to  have  an  initial  gaussian 
profile  falling  off  on  a  length  scale  of  500  meters. 

. 

nH2o+  =  nH]0+(z  ~  0)  exp 

where  |r,j  =  50  meters  and  Ri  -  500  meters.  This  simulates  a  large  scale  ionic  water  injection 
The  initial  water  ion  profile  for  a  typical  simulation  is  shown  in  Figure  3  where  the  central  density 
for  this  case  waa  nHj0+(x  =  0)  =  106  cm"5.  In  this  figure  the  x-direction  corresponds  to  the 
negative  wind  velocity  direction.  That  is  the  neutral  wind  is  blowing  from  positive  x  to  negative  x 
in  this  figure.  The  negative  y  direction  is  the  direction  of  the  motional  electric  field. 

In  Figures  4-6  we  show  the  initial  distribution  of  potential  for  a  density  profile  as  given  in  Figure 
3  and  for  three  initial  ion  water  densities.  In  Figure  4  the  initial  water  density  was  n//aC,  +  (f  = 
0)  =  104  cm"5.  In  this  figure  we  see  the  shielding  that  arises  when  the  total  ion  density  exceeds 
the  ambient  oxygen  ion  density  and  the  cloud  behaves  as  a  classical  dielectric.  In  Figure  5  the 
central  water  density  was  taken  to  be  nffj0+(x  =  0)  =  106  cm"3.  Here  we  see  exactly  the  opposite 
tendency  from  Figure  4.  The  electric  field  inside  the  cloud  increases  rather  than  decreases.  This 
occurs  because  the  density  gradient  terms  in  the  current  balance  equation  are  large  enough  to 
overcome  the  natural  polarization  which  occurs  due  to  the  imposed  motional  electric  field.  The 
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electric  field  increases  since  a  large  mobility  driven  electric  flux  is  necessary  to  balance  the  density 
gradient  driven  ion  flux.  This  tendency  is  confirmed  in  Figure  6  where  we  give  the  potential 
distribution  for  i  =  0)  =  108  cm'5.  The  potential  contours  here  can  be  well  modeled  as 

due  to  a  uniform  motional  field  and  a  positive  charge  on  the  positive  y  side  of  the  cloud  and  a 
negative  charge  on  the  negative  y  side  of  the  cloud.  This  suggests  that  the  picture  of  the  plasma 
cloud  around  large  objects  in  the  ionosphere  as  being  a  dielectric  shield[Katz  et  al.,  1984]  is  too 
simplisitic  a  model.  The  important  effect  of  the  plasma  density  gradients  must  also  be  taken  into 
account. 

A  typical  simulation  is  shown  in  Figures  7-13  for  the  initial  water  ion  puff  given  in  Figure  3  and 
with  an  initial  density  of  nH7Q+(x  =  0)  =  10®  cm-3.  Since  an  enhanced  density  of  this  magnitude 
would  have  to  be  caused  by  some  means  such  as  a  plasma  source  which  would  probably  emit  warm 
electrons  we  took  the  recombination  rate  for  the  water  as  kreeom  ~  3  x  10-8  cm8/secs[Murad  and 
Lai,  1986b]  which  corresponds  to  an  electron  energy  of  1  eV.  The  neutral  water  is  assumed  to  have 
the  form  given  in  Eq.  (1).  In  Figure  7  we  show  the  density  contours  at  0.1  seconds  after  the  start 
of  the  simulation.  This  corresponds  to  16  gyroperiods  which  is  enough  time  so  that  a  substantial 
E  x  B  drift  can  occur.  The  water  ion  cloud  is  seen  to  be  drifting  backwards  and  to  be  undergoing 
some  distortion  on  the  backside  of  the  cloud.  The  central  region  of  the  cloud  is  seen  to  be  breaking 
into  two  regions  with  very  steep  gradients  on  the  backside  of  the  central  region.  The  fact  that 
most  of  the  distortion  is  concentrated  on  the  backside  of  the  cloud  can  be  understood  by  noting 
that  the  high  density  core  of  the  cloud  will  drift  more  slowly  than  the  rear  fringes  of  the  cloud 
so  that  the  back  edge  of  the  cloud  will  collide  with  the  central  core  so  leading  to  the  observed 
distortion.  Also  around  the  source  ionic  water  is  still  being  formed  due  to  charge  exchange  and 
then  immediately  drifting  backward  so  that  a  very  steep  density  gradient  in  the  ionic  water  is 
observed  near  the  source.  The  distortion  on  the  backside  of  the  cloud  may  be  the  development 
of  the  E  x  B  instability.  In  barium  cloud  studies  where  the  equationsjPerkins  et  al. ,  1973]  for 
the  barium  cloud  are  very  similar  to  the  ones  used  here  the  cloud  was  observed  to  break  up  into 
filaments  both  experimentally  and  numerically  [Zabusky  et  al,  19731.  For  barium  clouds  it  is  believed 
that  the  filaments  may  be  associated  with  the  formation  of  E  x  B  instabilitiesjDrake  and  Huba, 


1986].  This  leads  us  to  suggest  that  the  filaments  may  be  due  to  a  similar  reason.  This  important 
topic  of  possible  instabilities  is  the  subject  of  ongoing  research  and  will  be  reported  in  a  future 
publication. This  may  be  another  reason  why  no  steady  state  may  exist.  This  also  suggests  that 
electrostatic  noise  may  be  inextricably  linked  with  the  large  water  densities  observed  around  large 
space  objects.  In  Figure  8  we  show  the  associated  oxygen  density  contours  at  the  same  time.  This 
region  delimits  the  area  where  new  ionic  water  is  being  created  from  the  area  where  ionic  water 
exists  because  it  is  came  from  the  initial  conditions.  Comparision  of  this  figure  with  the  last  figure 
suggests  that  part  of  the  breakup  of  the  central  core  may  be  due  to  the  creation  of  ionic  water 
there  in  an  asymmetric  manner.  In  this  figure  we  also  see  the  very  sharp  density  gradients  in  the 
oxygen  density  along  the  edges  of  the  oxygen  hole.  In  0.1  seconds  the  region  of  oxygen  depletion 
extends  back  behind  the  structure  a  distance  of  approximately  0.7  km.  This  suggests  that  the 
center  of  the  water  cloud  is  moving  with  a  drift  velocity  of  approximately  7  km  per  second.  This 
is  less  than  the  velocity  an  individual  ion  would  drift  backwards  with  if  it  saw  only  the  motional 
electric  field  and  occurs  because  the  ions  in  the  cloud  are  being  affected  by  the  self-consistent 
electric  field  and  not  the  motional  electric  field.  In  Figure  9  we  show  the  potential  contours  at  0.1 
seconds.  The  fact  that  the  cloud  has  moved  relative  to  the  structure  is  clearly  seen  in  this  figure 
since  the  region  where  the  potential  contours  is  distorted  has  moved  backwards.  For  Figure  10-13 
we  show  the  time  development  of  the  water  ion  cloud  at  0.15,0.2,0.25  and  0  3  seconds.  The  cloud 
continues  to  develop  a  tadpole-like  shape  in  the  direction  of  cloud  motion  and  the  central  breakup 
and  gradients  become  more  pronounced.  The  instability  on  the  backside  of  the  cloud  has  entered 
a  nonlinear  state  where  distinct  fingers  have  formed  on  the  backside  of  the  cloud.  This  has  also 
been  seen  in  barium  cloud  releases  in  the  ionosphere.  The  cloud  distribution  over  this  period  of 
time  suggests  that  measurements  of  the  cloud  density  from  the  structure  will  see  a  density  that  is 
both  time  dependent  and  spatially  anisotropic.  This  suggests  that  density  measurements  from  the 
Space  Shuttle  must  be  interpreted  with  care. 

The  time  that  ions  spend  in  the  vicinity  of  the  structure  is  critical  to  determining  the  ranee 
of  chemical  reactions  that  they  can  undergo.  This  question  is  of  substantial  interest  since  there  is 
evidence  that  chemical  reactions  can  occur  over  surfaces  in  space  Green  et  al  ,1985  In  Table  1 
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Table  1:  Residence  time  against  central  ion  water  density 

we  give  an  estimate  of  the  residence  time  of  a  water  ion  in  the  vicinity  of  the  Space  Shuttle  as  a 
function  of  the  initial  central  water  ion  density.  The  residence  time(r)  was  defined  as  r  =  L/vq 
where  L  is  the  length  of  the  Shuttle  taken  as  50  meters  and  vq  is  the  cloud  drift  velocity  determined 
from  the  simulation.  We  see  that  the  residence  time  of  the  ions  is  in  the  millisecond  to  tens  of 
milliseconds  range  and  futhermore  the  higher  density  cases  lose  ions  substantially  faster  than  the 
low  density  cases.  This  suggests  that  measurement  of  the  decay  rate  of  a  puff  of  plasma  will  not 
give  a  unique  answer  for  the  ion  residence  time[Sasaki  et  al.  1985].  Rather  the  ion  residence  time 
is  a  function  of  the  density  as  we  would  expect  for  a  loss  process  which  is  not  linear.  For  the  low 
density  case  the  residence  time  is  consistent  with  measurements  from  the  Shuttle;Caledonia  et  al.. 
1987],  In  Caledonia  et  al[l987|  it  was  shown  that  a  signal  of  n//jC)+  /n0  +  of  0.1  was  consistent  with 
a  loss  time  for  the  water  ions  of  40  milliseconds.  We  also  note  that  these  times  are  consistent  with 
results  from  Spacelab  2  where  enhanced  ionic  densities  were  not  seen  and  substantial  shielding  of 
the  motional  potential  was  not  observedlJ.  Raitt,  private  communication,  19861. 

5  Conclusions 

We  have  formulated  a  simple  model  for  the  plasma  cloud  around  a  large  structure  in  the 
ionosphere.  One  such  structure  is  the  Space  Shuttle  and  another  may  be  a  system  like  the  Space 
Station.  The  model  has  been  applied  to  an  ionic  water  puff  and  the  dynamic  behavior  of  the  cloud 
followed. 

A  number  of  important  conclusions  can  be  drawn  First  is  that  the  cloud  possesses  no  steady 
state  if  the  electron  density  is  enhanced  above  the  ambient  The  second  is  that  the  effect  of  density 


gradients  in  the  cloud  may  actually  enhance  the  electric  field  that  the  cloud  sees  so  that  the  picture 
of  the  cloud  always  acting  as  a  shield  against  the  motional  potential  is  too  simplistic  The  third  is 
the  cloud  may  be  subject  to  electrostatic  instabilities  which  will  grow  and  saturate  with  distortion  of 
the  cloud  structure.  This  will  mean  that  electrostatic  noise  and  turbulence  may  always  accompany 
these  clouds.  Finally  we  conclude  that  the  cloud  structure  will  be  highly  anisotropic  which  suggests 
that  measurements  of  the  cloud  density  from  the  structure  must  be  treated  with  care 

Future  work  will  concentrate  on  elucidation  of  the  electrostatic  instabilites  associated  with  the 
cloud  and  on  resolution  of  the  ion  residence  times  with  the  observed  chemistry  around  the  Shuttle 
These  conclusions  indicate  that  sensors  on  board  systems  like  the  Space  Shuttle  or  a  Space 
Station  will  have  to  cope  with  a  background  which  is  spatially  anisotropic,  temporally  varying  an  i 
also  electrostatically  noisy.  These  issues  will  also  be  studied  in  future  work. 
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We  substitute  Eqs.  (25)  and  (26)  into  Eq.  (20)  and  the  equivalent  equation  for  the  water  ions  and 
use  the  definitions  of  the  potentials  in  Eqs.  (23)  and  (24)  to  obtain 

-  V;  •  Vn0+  =  -^(V^0+  x  6)  •  Vn0+  (27) 

and 

-  Vc  ■  VnHj0+  =  -g(V^/f3o+  x  6)  •  Vnffj0+.  (28) 

These  equations  have  the  general  solution 


V0 o +  —  —H  + - 

no+ 


V'pHiO*  -  + - — ^*#30  + 

n//,o+ 


where  H  and  G  are  constant  vectors  related  to  Vc  and  and  are  arbitrary  functions 

of  the  water  ion  density  and  oxygen  ion  density  respectively.  If  we  substitute  these  expressions  for 
the  potentials  in  Eq  (22)  we  obtain 

V2('fr0+  4-  'I ~  M  ■  Vrj(j+  +  S  ■  VnHiO*  (^1) 

where  the  vectors  M  and  S  are  independent  of  the  densities.  We  integrate  Eq.  (31)over  a  volume 
V  bounded  by  a  surface  of  constant  electron  density  to  obtain 


/^V2(*0.  -  -0 


Now  if  we  multiply  Eq  (3!)  bj  integrate  over  the  same  volume  and  use  Eq  (32) 

and  Green’s  theorem  we  obtain 


This  is  only  possible  if 


(V*,,.  -  (1 


(V*„.  ‘  v*WjI).)  0 


Vji  *  M  J,  JiV-  ^  a  Vm  !L.\iiVk  Vm  'A 


\  .  ■  N  •  * 


This  indicates  that  the  electron  density  equation  must  be  of  the  form 

1?c.Vne  =  0. 

The  only  nonconstant  solution  to  this  equation  is  for  the  electon  density  to  vary  in  only  one  directiun 
in  space  which  is  physically  unreasonable.  Therefore  we  can  conclude  that  if  the  electron  density 
has  closed  contours  then  it  has  no  steady  state  and  if  there  is  a  steady  state  then  the  electron 
density  must  be  the  ambient  density  everywhere. 
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